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Abstract 

We investigate the scalar perturbation in the Lee- Wick bouncing universe driven by an 
ordinary scalar field plus a ghost field. We consider only a symmetric evolution of the uni- 
verse and the scalar fields about the bouncing point. The gauge invariant Sasaki-Mukhanov 
variable is numerically solved in the spatially flat gauge. We find a new form of the initial 
perturbation growing during the contracting phase. After the bouncing, this growing mode 
stabilizes to a constant mode which is responsible for the late-time power spectrum. 
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1 Introduction 

Lee and Wick (LW) constructed quantum electrodynamics with the higher derivative (HD) propa- 
gator for a photon [T] . The HD term in the denominator of the propagator improves the ultraviolet 
convergence of the Feynman diagram since the propagator falls off more quickly as momentum 
grows. A minimal set of the HD term leads to an additional physical pole in the propagator, 
corresponding to a massive LW-photon with the wrong-sign residue. It seems that the wrong sign 
gives rise to the instability or the violation of unitarity in the theory. Lee and Wick proposed a 
deformation of integration contours in the Feynman diagram so that the theory can be free from 
the instability and the violation of unitarity. 

Recently Lee- Wick Standard Model (LWSM) has been proposed as a candidate to solve the 
hierarchy problem in the Standard Model [2]. Every field in the LWSM has a higher derivative 
kinetic term. Using the auxiliary field (AF) method, the HD term can be converted into the 
degree of freedom of a massive field with wrong-sign kinetic term, referred as the LW partner. 
Since the LW partner can decay into ordinary fields, the wrong sign does not cause the violation 
of unitarity at macroscopic scales [HElll]. Soon after, the = 3 HD theory was constructed by 
mapping among the HD Lagrangian, the AF Lagrangian, and the LW form in Ref. [5], and the 
generalized formalism for the scalar field was constructed in Ref. [B]. 
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Very recently, Cai et. al. applied the N = 2 scalar-field Lee- Wick model to cosmology, 
and investigated the evolution of the universe [7j. (See also Refs. [8].) They found that the 
universe naturally bounces from the contracting phase to the expanding phase owing to the energy- 
condition violating Lee- Wick partner field. The merit of the bouncing universe in general is that 
cosmological problems such as the flatness and the horizon problems are naturally solved [7]: 
during the contracting phase, the deviation of Qk from decreases (solution to the flatness 
problem), and the universe starts well inside the horizon (solution to the horizon problem). The 
remaining is whether or not, the bouncing universe can provide proper density perturbations for 
structure formation. What the authors found in the Lee- Wick bouncing model are that the density 
perturbation at the bouncing point is nonsingular and the power spectrum can be scale-invariant. 
See also Refs. [9] for the recent works about nonsingular bouncing universe models. 

In this work, we investigate the scalar perturbation and study precisely its properties by 
inspecting the comoving curvature. We consider a symmetric background universe about the 
bouncing point and numerically solve the scalar perturbation. We find a new type of initial 
vacuum solution, and discuss its growth and resulting late-time power spectrum. 

The paper is organized as following. In Sec. 2, the evolution of the background universe is 
calculated. In Sec. 3, the scalar perturbation theory is briefly introduced. In Sec. 4, the gauge 
invariant Sasaki-Mukhanov variable is calculated. In Sec. 5, the comoving curvature is evaluated, 
the initial perturbation is obtained, and the late-time power spectrum is discussed. In Sec. 6, we 
conclude. 



2 Evolution of background universe 

The general higher derivative Lagrangian of a self-interacting real scalar field (p is given by 

1 ^ 1 

n=l 

where □ = -^^d^ {^—QQ^^dJ) , denotes the number of physical poles in the propagator of </), 
ttn is a constant coefficient with mass dimension [a„] = 2 — 2n, and V^(0) represents an interac- 
tion potential of the scalar field. Introducing auxiliary fields we can show the Lagrangian (12.11) 
is equivalent to that of scalar fields with quadratic kinetic terms, referred as the Lee- Wick 
Lagrangian, up to the quantum level [6], 

1 ^ 

4w = 2 E(-l)"^Vn(n - mD^r. - V(^i, ■ ■ ■ , ^^), {ml<ml---<ml), (2.2) 

n=l 

where V represents the interaction potential originated from V in (12. ip . 
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2.1 Symmetric bouncing universe 



We consider the N = 2 Lee- Wick model with vanishing interaction potential, which consists of one 
ordinary massive scalar field (fi and one ghostlike Lee- Wick partner field (f2- The Einstein-Hilbert 
action with these matter fields is 



^1 



LW 



N 



LievrG 



n=l 



(2.3) 



The corresponding energy- momentum tensor is given by 

N 

'2 



(2.4) 



n=l 



Nonvanishing components of the energy-momentum tensor for the time-dependent scalar fields 
are 

AT 



n=l 
AT 



n=l 



1.0 1 22 



(2.5) 



With an isotropic metric ansatz describing the background universe, 

ds'^ = -df + a{tfdx'dx\ 
the Einstein's equation and the scalar field equation read 

N 



(2.6) 



ij2 = ^E(-i; 

n=l 



N 



H = -47rGE(-l)"+V 



n=l 



(2.7) 

(2.8) 
(2.9) 



Since the Lee- Wick partner field is a ghost that violates the null energy condition, the background 
evolution can have a bouncing from contracting to expanding. Let us set the bouncing moment 
at t = 0. In addition, we shall consider a symmetric background about t = 0. In order to have a 
symmetric bouncing universe, we need to set conditions at t = 0. The choice of those conditions 
is not unique, and we adopt one as following: 
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1. One should have a(0) = 0, i.e., H{0) = 0, as a necessary condition for a bouncing universe. 
This makes the left-hand side of (12. 7p vanish. One can arrange the right-hand side in several 
ways to make it vanish. We adopt nonvanishing values of v?n(0), then V9„(0) should vanish 
for (pn(t) to be symmetric, i.e., (pnitys are even functions. 

2. As v?n(0) = 0, H{0) = from ^M)- This implies a(0) = 0. 

3. Differentiating fl2.8p . it is manifest that a (0) = 0, i.e., H{0) = 0. 

4. The conditions 1-3 provide a bouncing universe. However, for a bouncing from contracting 
to expanding, another condition is needed, 'a'(O) > 0, i.e., H{0) > 0. In order to satisfy this 
condition, we shall see later that one must have |mi| < \m2\. 

As a whole, a{t) and (pn{t) are even functions, and H{t) is an odd function. 
2.2 Numerical solutions 

We numerically solve the field equations (12. 8p and (12.91) . We integrate the field equations from 
t = to the t > region. What is good to do this is that it is numerically convergent and stable. 
Then the solutions for the t < region are obtained simply by symmetry. If one integrated from 
some moment at t < 0, the numerical instability would grow up, and it should be very difficult 
to overcome the bouncing point. In addition, even if one could obtain smooth regular solutions 
connected to the t > region, it would never be strictly symmetric about the bouncing point. 

Applying the symmetric bouncing conditions at t = discussed in the previous subsection, 
once m„'s are fixed, the only condition remaining for numerical calculation is v^i(O), or (p2{0). The 
numerical results for the background fields are plotted in Fig. 1. For numerical calculations in 
this paper, we set m2 = 2mi = 10~^mp = lOM^ and (fi{0) = {fn2/mi)(f2{0) = f^p = lO^MrE 

At the bouncing point t = 0, the scalar fields start to decrease from their maximum values, 
and then undergo damped oscillations about ipn = 0. In the t < region, the evolution is exactly 
opposite; the amplitudes of oscillations grow and the scalar fields reach their maxima at t = 0. 
The amplitude of ipi is much larger than that of ip2 in our set-up. 

While the scalar fields roll down from their maxima for t > 0, the universe undergoes accelerat- 
ing expansion, a > 0. Afterwards, the scalar fields start to oscillate and play the role of pressureless 
dust which drives the universe to settle down to the matter-dominated expansion, a ~ 

t2/3. When 

t < 0, the universe contracts initially as a ~ (— t)^/'^, and then undergoes accelerating contraction 
till the bounce. 

^In this paper, we introduce a mass unit AI^ = 10~^mp, where rup is the Planck mass, and aU the physical 
quantities are expressed in this unit. 
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Figure 1: Plot of the background fields which are symmetric about t = for m2 = 2mi = 10 ^rup = 
WMr, ifi{0) = (m2/mi)v?2(0) = rup = lO^M^, (^i(O) = (^2(0) = 0, H{0) = 0, and a(0) = 1. (a) Scalar 



fields ifi and (p2 exhibit damped oscillations. The amplitude of (pi is much larger, (b) The scale factor 
a shows the accelerating contraction/expansion near t = 0, and contraction/expansion at large 
(c) The horizon scale diverges near t = and at large 
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The horizon scale \H~^\ diverges at t = and t — ?■ ±00. A moderate comoving scale initially 
starts well inside the horizon at large \t\ (t < 0), and crosses the horizon four times to enter the 
current universe. 

During the matter-dominated period controlled mainly by oscillating ipi at large \t\, the scale 
factor is approximated by a ^ at^^^ where a is a constant depending on the values of m„ and 
V9„(0). Then field equations fl2.7p -( l2l9|) allow the following asymptotic solutions, 

cos(mit + ai) 
V ottG m\t 

m-^-^"'"^^""''\^'"'\ (2.11) 

where ai is a phase factor depending on the values of m„ and ifniO)- 



3 Scalar Perturbation and Sasaki- Mukhanov Variable Qn 



In this section, we briefly review the scalar perturbation with multiple matter scalar fields, based 
on Refs. [T0] - [T2] . and apply for the case of N Lee- Wick scalar fields. The metric perturbation is 
decomposed into four scalar, four vector, and two tensor components. For the scalar perturbation 
with matter scalar fields, the total number of perturbation modes is 4 + A^. Imposing two 
constraints originating from (00)- and (0«)-components of the Einstein's equation and fixing two 
gauge degrees of freedom from the (4 + A^)-modes, we have only A^ dynamical degrees of freedom 
in the system. 

Let us consider the linear scalar perturbation of the Friedmann-Robertson- Walker metric, 

ds"^ = -(1 + 2A)dt'^ + 2adiBdx'dt + [(1 - 2ip)5ij + 2didjE] dx'dx\ (3.12) 

where A, B, ip, and E are four scalar modes. Then the (00)- and (02)-components of the perturbed 
Einstein's equation give the relations, respectively. 



SHiilj + HA) + — iP + Hia'E - aB) = -AiiGSp 
^ + HA = -AirGSq, 



(3.13) 



where k stands for the comoving wave number, and 

N 



-1 



in+l 



n=l 



N 



(3.14) 



n=l 
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and 5(pn represents the perturbed mode of the background scalar field The field equation for 
the matter fields is given by 



a^E - aB 



(3.15) 



After fixing two gauge modes, we can completely determine the time evolution of the scalar 
perturbation using the equations fl3.13p and fl3.15p . In other words, the resulting behaviors of the 
scalar perturbation satisfy the remaining (i,j)-components of the perturbed Einstein's equation. 

In order to fix two gauge degrees of freedom, we introduce several gauge invariant quantities 
in terms of the scalar modes in the metric (I3.12p . 



^ = A + aH{B - aE) + a[B - aE)' 
^ = ij-aH{B -aE). 



(3.16) 



We also introduce a gauge invariant definition, which is refereed as Sasaki-Mukhanov variables |13j . 

(3.17) 

In this paper we choose the spatially fiat gauge. 



^ = 0, E = 0. (3.18) 
Then we can replace gauge dependent variables in (13.131) and (I3.15P with the gauge invariant ones, 

S'-Pn = Qn, 



a^E - aB 



H' 



(3.19) 



Using these relations we can rewrite the constraints in f l3.13p as 



^ = -A7rG{5p - 3H5q) 



N 



n=l 



-47rG^(-l)"+l \MQn - ^nA) - (pnQn 



N 



A 



-- = -Jf 2^(-l)^VnQn. 

71=1 



(3.20) 



Plugging fl3.19p and fl3.20p into fl3.15p . we express the right-hand side of f l3.15p with the background 
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variables and Qn, 



RHS = -2m>„A + Lpn[A + 



N 



1=1 



dt \H 



8ttG 



1=1 ^ ^ 



(3.21) 



where we used the background evolutions (12.81) and (12.91) . We finally obtain N equations for 



Qn + 3i/Q„ + — Q„ + mlQn - ^ ^(-i)'+i_ ( ]Qi = 0. 

1=1 



dt\H 



With the Sasaki- Mukhanov variables Qn, the comoving curvature TZ is defined as 



■R = iIj 



H 



P + P 



5q = H X 



ELi(-i)"^-^V 



(3.22) 



(3.23) 



4 Evaluation of Qn 



In this section, we solve numerically (I3.22p for = 2 to obtain Qn{t). The background fields ^n{t) 
and H{t) were already solved numerically as in Sec. 2. In order to solve Qn{t) in this background, 
first we perform the series expansion about the bouncing point t = 0. The background fields can 
be expanded as 



^nit) = PnO + Pnlt + Pn2t^ + PnS^^ + 

Hit) = ho + hit + h2t'^ + ht^ + ■■■ 



(4.24) 
(4.25) 



Plugging these expansions in the background field equations (12. 8p and (12.91) . we can obtain the 
coefficients Pni^s and hi's. As we mentioned in Sec. 2, ip^s are even functions and H is an odd 
function. Therefore, Pniii = odd) = and hi{i = even) = 0. In addition, applying the symmetric 
and bouncing conditions at t = 0, all the other coefficients are expressed in terms of Pno- [Note 
that pno = V5n(0). For the N = 2 case, only one of Pno is free since m^y9^(0) = m^v^llO) from (12.71) .] 
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The coefficients are then determined as follows; 

,2 



Pn2 
PnA 



24 



Pn6 = - 



30 



N 4 1 ■ 



PnO, 



(4.26) 



AttG 



N 



3 

47rG 



j2i-^r^'<pio^ 



n=l 
N 



15 



^(-l)"+^m^p: 



2 



n=l 



(4.27) 



The corresponding behavior of the scale factor near the bouncing point is given by 

a{t) = 1 + ]hst^ + Iht'^ + ■■■ 
4 



(4.28) 



with the normalization a(0) = 1. As it was pointed out in Sec. 2.1, for the bouncing from 
contraction to expansion, one should have 'a'(O) > 0. This is satisfied when > which implies 
\m2\ > |mi| using (14.271) . 

Now we expand Q„ in series as 



Qn(t) = f{qnO + Qnlt + g„2t^ + Qn^f + qnAt"^ H )• 



(4.29) 



Since the comoving curvature IZ in (I3.23P can be regular even when is singular, we introduced 
factor in the front (and it can be s < 0). Plugging the expansions (I4.26p . fl4.27p . and (I4.29P 
into the Q-equation (I3.22p . we obtain two linearly independent solutions for Q„. The one is an 
even function with s = —2 and g„j(i = odd) = 0, and the other is an odd function with s = and 
qni{i = even) = 0; 



Q--(t) = t~\qno + g„2t' + qn,t' + ■ 
QTit) = Quit + qnst^ + Qn^t' + ■■■ . 



(4.30) 
(4.31) 
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In addition, we are also provided with the relations for for the even case, 

920 = — 9io, (4.32) 
mi 

mi {mj - ml){5k^ + ml + ml) 

922 = — 912 910, 4.33 

m2 3Umim2 

and for the odd case, 

921 = — 911, (4.34) 
m2 

mi Dmim2 

All the other g„j's can also be determined by the background parameters (mi, m2) and the 
fluctuation parameters (/c, gio, 912, 9ii, 913)- For a given comoving wave number k, qio and gi2 
(gii and ^13) are free parameters for the even (odd) case. These free parameters are the shooting 
parameters in numerical calculations for Q„. Since the Q-equation is linear, constant xQn is also 
a solution. One of the shooting parameters plays a role of constant multiplication, and the other 
plays the role of real shooting parameter. From numerical calculations, we observe that the shape 
of the solutions does not depend on the value of the second parameter, which means that the 
solutions are attractive. 

The numerical solutions of Qn for the even and the odd cases are plotted in Fig. 2. As 
the background scalar fields ipnS enter the oscillating regime, the universe settles to the matter- 
dominated expansion. Then the perturbation Qn oscillates about zero apparently with a constant 
amplitude in the same frequency with the background fields. The amplitude of Qi is much larger 
than that of Q2, so the perturbation properties such as the comoving curvature TZ are mainly 
determined by Qi- 

For future purpose, let us analyze the late-time behavior of the perturbation. Since Qi is 
dominant, neglecting Q2 in (13.221) . the perturbation equation for Qi is approximated by 

Qi + SiJQi + + ml^ Qi - ^1 (^vi?) Qi ^ 0. (4.36) 

Plugging the asymptotic solutions for (fi and H presented in (I2.10p and (12. lip , and a{t) ~ at^/^, 
and selecting the dominant terms at late times {mit ^ 1), the above equation becomes 



2 ■ 

Q1 + -Q1 + 



c,2-^ + - ^ sin(2mit + 2ai) 



Qi ^ 0. (4.37) 



Among five terms in the equation, the first term and the fourth (mass) term are the most dominant 
and comparable in magnitude. These terms give a constant oscillation. The next dominant terms 
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Figure 2: Plot of the Sasaki-Mukhanov variable Qn for the even case {qio = 10~^^, qi2 = 10) and for the 
odd case {qn = 10^, gis = 10^). Qn settles down to an oscillation with an almost constant amplitude. 
Qi has a much larger amplitude. The t < region is evenly and oddly symmetric. For the even case, Qn 
diverges at t = 0, but is not visible in the figure. 

are the second (friction) term and the last (sinusoidal) term which are comparable in magnitude. 
These two terms add damped oscillation behavior in Qi. The third (k-) term is comparable 
with these two terms in magnitude during the intermediate period for large k. In this period, the 
damped oscillation in Qi behaves as t~^^^ due to the /c-term. The numerical results shown in Fig. 2 
contain the constant and the damped oscillations in principle, although the damped oscillation 
is not visible because it is tiny compared with the constant oscillation. These two oscillating 
patterns play an important role in the comoving curvature TZ as we shall see in the next section. 

5 Initial Perturbation and Power Spectrum 

5.1 Comoving curvature IZ 

For N = 2, the comoving curvature TZ in fl3.23p becomes 

n = {ipiQi - 02Q2) = fiQi - /2Q2. (5.38) 

The factor /„ is a function of the background fields, and is an even function. Therefore, the 
evenness/oddness of TZ is the same with that of Qn- The factor function becomes a periodic 
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Figure 3: Plot of the factor function fi and fiQi in TZ. (Only the odd case is plotted, but the even 
case is similar.) The shape of /2Q2 is also similar, but the overall level is lower and the period is half. 



function at late times. (See Fig. 3.) Whenever ipl = /„ diverges. However, this divergence 
does not make the evaluation of TZ unphysical in the whole time (THUS] . One may ignore simply 
the diverging period only while keeping the rest as physical. 

As we observed in the previous section, Qn has a periodic oscillating behavior. The frequency 
is the same with that of but the phase is a little bit off. Then the shape of fnQn becomes 
as shown in Fig. 3. Considering one period, except for the vicinity of the divergence, fnQn stays 
almost constant. In the successive periods, this pattern is maintained, and the value of fnQn does 
not seem to change due to the apparently constant amplitude of Qn- However, as it was discussed 
in the previous section, there is a damped oscillation in Q„ in addition to the constant oscillation. 
The effect of this damped oscillation will be implied in the pattern of fnQn- Therefore, as periods 
proceed, the value of fnQn will decrease slightly. 

The behavior of fnQn is embodied in the comoving curvature TZ from f l5.38p . As it is shown in 
Fig. 4, TZ stays constant except the vicinity of the divergent points, while the effect of the damped 
oscillation is not very visible. 

What is remaining is how to extract the damped oscillation from our numerical results of 
TZ- In order to see this, we plot log^Q |7?.| vs. log^o^ in Fig. 5. The slope is ~ —1/3, which 
indicates that the decaying hidden in constant TZ behaves as T^^'^^aymg ^ t^^/^. This decaying 
behavior results from the damped oscillation in Qi controlled by the fc-term in fl4.37p discussed 
in the previous section. In the numerical domain shown in the plot, the effect of the fc-term is 
considerable. If the numerical calculation could go further in time, one would be able to see the 
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logio \n 



Even 



Odd 



(a) (b) 

Figure 4: Plot of the comoving curvature TZ. The spikes indicate divergences. Except for the divergent 
region, the curvature approaches an apparently constant value soon. The constant value can be adjusted 
by varying the values of shooting parameters. The decaying in IZ from the damped oscillation in Qi is 
inherent, but not visible. 





logio t 



Figure 5: Plot of log^^g 1^1 vs. log^gt. The slope indicates that IZ contains the decaying mode. (We 
plotted only for the odd case, but the even case is similar.) From the slope, the decaying is proportional 
to t~^/^ . The decaying behavior is not achieved within the numerical domain. 
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decaying behavior transit from t to t ^. Unfortunately, however, the numerical integration 
stops before the transition due to numerical error accumulated. 

5.2 Initial Perturbation 

Since the comoving curvature TZ has been obtained, let us discuss the initial perturbation consid- 
ering the modes of TZ. In the previous section, we obtained two linearly independent modes of TZ, 
the even {E-) and the odd (0-) modes. In general during the expanding phase (t > 0), the curva- 
ture consists of the constant (C-) and the decaying {D-) modes which are linearly independent. 
The C- and Z^-modes can be obtained from a linear combination of the E- and 0-modes. These 
C- and D-modes are those discussed in the previous subsection. 

During the contracting phase (t < 0), the decaying mode becomes a growing (G-) mode due 
to symmetry. The G-mode provides the initial perturbation in the bouncing universe. The initial 
perturbation is produced in the contracting phase deep inside the horizon, and then grows. The 
perturbation produced in this way is frozen after it crosses the horizon in the expanding phase 
(t > 0), and should produce the observed value in the power spectrum. 

In order to discuss the initial quantum fluctuation and the vacuum state, let us introduce the 
conformal time r] and a new field v for canonical quantization by 



We consider only the contribution of Qi which is dominant at late times as discussed in Sec. Si 
The functional behavior of v at early times [rj <^ 0) can be read from that at late times (r/ ^ 0), 
due to the symmetric property of our setting. The action for the new field is written by 



dt = adr], 



V = aQi. 



(5.39) 




(5.40) 



where ' = 



d/drj, z = a(p/H. The general solution of the equation of motion is given by 

~ / ^^ f d^k \ , , ^ ^, , N + 1 ih.:^ 




(5.41) 



where v{ri\ k) satisfies 




(5.42) 



and is normalized as 



vv 



* / 

— V V = I. 



(5.43) 



^We fix the bouncing point as rj{t = 0) = 0. 
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At early times, 77 <^ 0, the field equation (15.421) can be approximated by 



-V 



-7] sm 



27 



+ 2ai 



V{7]) ^ 0, 



(5.44) 



where we used a ~ at^^^ = a^rf I'd and t = o?rf jll . This equation can be obtained from (15. 42 p . 
or from a direct transformation of (I4.37p . The horizon crossing for a given wavelength occurs 
when 



A 



•ph 



H 



\kr]\ ~ 2. 



In the far subhorizon limit, \kr]\ ^ 2, the equation (I5.44p can be further approximated by 
The solution to this equation satisfying the normalization condition (I5.43P is given by 



(5.45) 



(5.46) 



v{ri) 




(1) / mia-" 3 



27 



V 



(2) / rnia'' 3 



27 



(5.47) 



where — l^il 



P = 1 and Hn'"^^ represent the Hankel functions. We adopt the Bunch-Davies 
vacuum for the production of the initial perturbation at 77 = 77^ <^ by taking only the positive 
energy mode, Ai = 0. Then the vacuum solution (15.470 at 77 ^ becomes 



v{r]) 



9 



2mia'^ ?7 



- exp 



—I 



rriia 3 vr 
27 ^ ~ 3 



(5.48) 



This solution represents the initial perturbation produced deep inside the horizon, and its ampli- 
tude grows in time. It can be transformed to the original variables at t ^ 0, 



^ , , V 1 . . . 
Qi{t) = - (X - exp [—tmit) 
a t 



(5.49) 



During the expanding phase (t > 0), this solution represents the damped oscillation analyzed 
in the previous section, which gives 7^'ie'='^y™g oc 

The functional dependence on 77 of the initially growing perturbation (I5.48P is somewhat differ- 
ent from the initial perturbation studied in inflation or in power-law expansion in the literature. In 
those, V ~ l/y/2k exp{—ikri). The difference originates from the oscillating behavior oi the back- 
ground fields in the Lee- Wick model, which provides the ry^-term in (I5.44p as the most dominant 
term. 
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5.3 Analysis of IZ modes and late-time power spectrum 

In this subsection, let us discuss the modes of TZ. The initial growing mode of 71 can be obtained 
from the results of the previous section. Since Qi has the dominant contribution to TZ, the G-mode 
of 7^ at t < is obtained from fl5^ . 

^growing _ J^Qgrowing ^ 1_ ^^ ^q^ 

This G-mode can be extracted from the E- and 0-modes by a linear combination. Each of the 
even and the odd mode contains the growing and the constant modes in itself. Then, the growing 
mode can be obtained as following; 

n{t < 0) = Cl7^"™° + C27^°<^'^ (5.51) 

^ l^-^cvcn— growing _|_ -^even— const j _j_ ^ j^-^odd— growing _j_ -^odd— constj 52) 

= 7^s™g, (5.54) 

where ci and C2 are fixed such that ci??.'^^'''^"™'^'^* + ^2?^°'^'^"™'^'^* = 0. Then the constant modes 
in TZ are canceled each other, and only the growing mode remains and plays the role of initial 
perturbation. 

Once one fixes Ci and C2 as above, the linear combination of 7^^™° and 7^°^^*^ is fixed. At t 3> 0, 
the comoving curvature 7?. as a result of this linear combination becomes 

7^(t > 0) = Cl7^'="'=" + C27^°'^'^ (5.55) 

^ l^-^cvcn— decaying _|_ -^cvcn— const j _|_ ^ j^-^odd— decaying _j_ -^odd— constj 56) 

= [cl7^°™"~'^'='^y^''^ + C27^°'^'^^'^°'''y'"s] + [cl7^'™''^™''''* + C27^°'^'^"™''''*] (5.57) 

^ ^decaying ^ ^const _ ^const^ gg^j 

As it can be seen in the schematic picture in Fig. 6, the constant modes in JZ^^^"^ and 7^°'^'^ do not 
cancel each other because the odd mode changes its signature at t > 0. Therefore, at late times 
during the expanding phase, the comoving curvature approaches a constant. This constant must 
reproduce the observed value of the power spectrum, 

Vn = l^l^l' ~ 10"'. (5.59) 

Although the schematic story is as above, it will be very difficult to extract the appropriate 
-^growing jjiode from numerical calculations. First of all, our numerical results could not reach 
the region of large enough \t\ where the decaying mode is realized. Even when the numerical 
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Figure 6: Schematic picture of the linear combination of T^'^™" and 7^° . At t ^ 0, the initial growing 
mode T^s™™™^ can be extracted from the combination. At t ^ 0, this combination provides the constant 
mode 7^^°^*^* which is responsible for the late-time power spectrum. 

technique is improved to reach such a large \t\ region, combining the even and the odd modes in 
such a way as to meet the both of the initial condition (15.481) and the late-time power spectrum 
fl5.59p simultaneously, is nontrivial. Since Qn has an oscillating motion, one needs to match the 
amplitude as well as the phase exactly by adjusting the initial shooting parameters. It is not clear 
whether or not the numerical error will admit such an exact matching. 

6 Conclusions 

In this paper, we investigated the scalar perturbation of the Lee- Wick bouncing universe driven 
by two scalar fields, one of which is an ordinary field and the other is a ghost field. The universe 
contracts initially and then expands, owing to the energy-condition violating ghost field. For the 
analytical purpose, we considered a symmetric evolution about the bouncing point t = 0; The 
scale factor a and the scalar fields ipnS are even functions. Near the bouncing point, the universe 
undergoes an accelerating contraction and expansion. If we set the ghost field is more massive 
than the ordinary field ipi, m2 > mi, the universe is mainly driven by the ordinary field at large \t\. 
The ordinary field (pi at large t exhibits a damped oscillation and plays the role of the pressureless 
dust driving the universe to a matter-dominated expansion. At t < 0, the evolution is opposite. 
We investigated the scalar perturbation in this background. For the scalar perturbation, we 
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adopted the spatially flat gauge and numerically solved the gauge invariant Sasaki-Mukhanov 
variable Qn- In calculating Qn, we performed a series expansion about t = 0, and found that 
there exist two linearly independent solutions. One is an even function and the other is an odd 
function. This kind of approach in calculating Qn about t = reduced difficulties in dealing 
with the bouncing point. (If one starts to solve numerically Qn from t < 0, it is very difficult to 
overcome the bouncing point.) The resulting evenness/oddnes of Qn coming from the symmetry 
of the background saved the labor for the half of calculation, and enabled us to analyze the 
corresponding physics in an easy way. 

From the numerical results of Qn-, we evaluated the comoving curvature TZ. From the linear 
combination of the even and the odd modes of 7^, we showed that the growing mode and the 
constant modes can be extracted in principle. The growing mode is responsible for the initial 
scalar perturbation produced deep inside the horizon during the contracting phase. We found 
that the initial perturbation in the Bunch-Davies vacuum has a form fl5.48p which is different from 
those known in the literature. This difference comes from the oscillating nature of the background 
fields. This growing mode is connected to the constant mode achieved after the perturbation 
crosses out the horizon during the expanding phase. This constant mode explains the late-time 
power spectrum. 

Due to the numerical difficulties, we could not see the decaying/growing behavior in TZ. 
With a better numerical technique, we wish to come back to achieve the complete decaying 
behavior in the future. We considered only for the N = 2 Lee- Wick model, but it will also be 
interesting to study for the N > 2 cases. 
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